Tristability in a non-equilibrium double-quantum-dot in Kondo regime 
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Electron tunneling through a non-equilibrium double 
quantum dot in the Kondo regime is studied. In the region of 
negative differential resistance, it is shown that this system 
possesses a complex response to the applied potential char- 
acterized by a tristable solution for the current. Increasing 
the applied potential or reducing the inter-dot coupling, the 
system goes through a transition from a coherent inter-dot 
regime to an incoherent one. The different nature of the so- 
lutions are characterized and it is shown that the effects of 
the asymmetry in the dot-lead coupling can be used to con- 
trol the region of multistability. The mean-field slave-boson 
formalism is used to obtain the solution of the problem. 



JQDs) has been 
□ The experimen- 



The Kondo effect in quantum-do 
extensively studied in the last years 
tal evidence has confirmed that many of the phenomena 
that characterize strongly correlated metals and insula- 
tors, as it is the case of the Kondo effect, are present 
in QDs. The QDs allow to study systematically the 
quantum-coherence many-body Kondo state, due to the 
possibility of continuous tuning of the relevant parame- 
ters governing the properties of this state, in equilibrium 
and non-equilibrium situations. 

The electron tunneling through double-quantum-dots ( 
DQD) in the K^^Qxeaki^lias received much attention 
in recent years. □ BO'Dffl ^'EJ'Ej In comparison with a sin- 
gle quantum-dot, the double-quantum-dot has a richer 
physics. For example in DQDs it is possible to study 
the competition between the inter-dot antiferromagnetic 
spin-spin correlation and the dot-conduction spin-spin 
correlation present in its ground state and also it can 
be investigated the inter-dot coherence effects that re- 
sults from the Coulomb interaction. The type of coupling 
between the QDs determines the character of the elec- 
tronic states and the transport properties of the artificial 
molecule. In the tunneling regime, the electronics states 
are extended across the entire system and can be con- 
structed from a coherent state based on the bonding or 
anti-bonding levels of the QDs. Some aspects of the non- 
equilibrium transport properties of a simplified DQD in 
the Kondo regime, constituted by two identical QDs has 
recently been studiedllil. When the inter-dot coupling 
is greater than the level broadening, there is a region of 
voltage where this system shows a bistable behavior char- 
acterized by two solutions for the current. This behavior 



resembles the J-V characteristic curve of a double bar- 
rier structure in the accumulatio n of charge regime and in 
doped superlatticesMiiS'tSrEyGdHowever, theoretical and 
experimental studies have confirmed that a double bar- 
rier structure possesses a more complex behavior than 
the one predicted by a bistability. The J-V curve of this 
system has in many cases a Z-shape, indicating the exis- 
tence of a tristable solution for the current O 



In this work we analyze the new physics derived from 
the complex behavior that characterize the electron tun- 
neling through a non-equilibrium DQD in the Kondo 
regime. Although this system is essentially different from 
a double barrier heterostructure, we find that, in certain 
region of the parameter space, the J — V characteristic 
curve of both are similar. However, the problem we ana- 
lyze has much richer physics. Increasing the asymmetry 
of the coupling of the dots with the leads or reducing the 
inter-dot coupling, the system goes through a transition 
from a coherent inter-dot regime to an incoherent one. 
The characteristic curve reflects this transition changing 
from a Z-shape in the incoherent region to a loop-shape 
when the dots behave coherently. The different nature of 
the solutions are characterized and it is shown that the 
effects of the asymmetry in the dot-lead coupling can be 
used to control the region of multistability. 

We adopt the two-fold degenerate Anderson Hamilto- 
nian in the limit U — > oo, that is diagonalized using the 
mean-field slave-boson formalism. Although this mean- 
field approximation has, in certain region of the param- 
eters that define the system, a non-physical solution it 
has been shown that the physics related to the-, nega- 
tive differential resistance is not associated to itO. The 
DQD is driven out-of-equilibrium by means of a dc bias 
voltage V at zero temperature. The double occupancy in 
each dot is forbidden and the inter-dot Coulomb inter- 
action neglected. We introduce the slave-boson opera- 
tor w a that creates an empty state and a fermion op- 
erator f a<7 that annihilates a single occupied state with 
spin a. To eliminate the possibility of double occupancy 
we impose the constraint Q a = J2 a flafaa + b\j3 a = 1. 
The annihilation operator of an electron in the dot a 
is c a<7 = b^facr. In the mean field approximation the 
bosonic operators b^ a and b a are replaced by their ex- 
pectation values, (b a ) = b a y2 = (fe^) = 6|,v2- Hence 
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the Hamiltonian of the double quantum-dot connected 
to leads plus constraints is written as, 



H — Hi ead + £an aa 

a=0,l,<7 



(i) 



a=0,l 



where £ _= s + X , si = ei + Ai, V L = V L b , V R = V R b\, 
t c = t c bobi,the X a are Lagrangian multipliers which guar- 
antee the constraint conditions on the Q a and the param- 
eters t c , Vl and V R are the inter-dot connection and left 
and right dot connections with the leads respectively. 

The first term on the right-hand side of Eq. 1 represents 
the electrons in the left and right leads; 



Hiead = E eiUia + 1 £ c \a c i"- ( 2 ) 



where the operator c\ a creates an electron in the site i 
with spin a, £j is the site energy and t is first-neighbor 
hopping in the leads. 

Here Hi ea d corresponds to the free-particle Hamilto- 
nian with cigenfunctions of the Bloch type, 



\k,a)=J2e ikia |J» 



(3) 



where \k,a) is the momentum eigenstate with spin a and 
\j) is a Wannier state localized at site j of spin a. The 
dispersions relation associated with these Bloch states 
reads 



£/c = 2tcos(ka) 



(4) 



The model has an energy band, extending from — 2t to 
+2t, with the first Brillouin zone expanding the k interval 
[— ir/a, n/a]. 

The stationary state of the complete Hamiltonian H 
with energy e k can be written as 



IV'fcer) 
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at \j, a) 



where the coefficients a\ a are given by, 
a k ja =< j, v\ipkcr > ■ 



(5) 



(6) 



We obtain the following eigenvalues equations for the 
Wannier amplitudes a^ a 

e k a) a = (j,a\H\ip ka ) 

e k a k ha = e 3 al a + t(a^_ ha + a k j+ha ) {j + -1,0, 1, 2), 

£ fe a -l(2) :C r = £-1(2)0-2(2), cr + ^(-R) a 0(l),cr + to -2(3),a-> 
£ fe a 0(l),(T = ^0(l) a 0(l),<r + ^L(fl) a -l(2),cr + *cai(0),ffi ( 7 ) 



where a k a is the amplitude of probability to find the 
electron in the site j and state k with spin a. 

The four parameters (bo, b\, Ao, Ai) are determined 
minimizing the expectation values of the Hamiltonian 
(7). They satisfy the set of four equations: 



5 o(d + ^E 

k.a 



a 0(l),o- 



Vi 



L{R) 



E Rc ( ab i(2)/S(i)J + 



E Re ( a l(0),<T a 0(l),<r) + ^0(l) b \(R) - °- 



(8) 



The sum over the wave vector k cover all the occupied 
states. The resulting equations are nonlinear because of 
the renormalization of the localized levels in the dots, the 
interdot coupling tunneling and the coupling tunneling 
between the QDs and the leads. 

The stationary scattering problem is solved by finding 
the eigenstates with positive k. In order to study the so- 
lutions of equations (4) and (5) we assume a plane wave 
incident from the left with an amplitude /, with a partial 
reflection amplitude R. The waveform at the right is a 
simple plane wave with intensity given by the transmis- 
sion amplitude T. Taking this to be the solution of the 
system, we can write, 



a) = Ie lka3 + Re- lka > ,j<0 



a' 



i k =Te lkRa \ j > 1 



(9) 



and those with negative k 



aj k = Ie- ikaj +Re ika \j > 1 

aj k = fe- lkLaj ,j < (10) 

These functions describe plane-wave electron ap- 
proaching the scattering potential from j = — oo and 
j = +oo, respectively, with transmission and reflections 
amplitudes T and R for k > 0, and T and R for k < 0. 

The solution of equations (7) can be obtained through 
an adequate iteration from right to left for k > 0, and 
from left to right for k < 0. For a given transmitted am- 
plitude, the associated reflected and incident amplitudes 
may be determined by matching the iterated function to 
the proper plane wave at the left lead for k > and 
at the right lead for k < 0. The transmission amplitude 
t(k) = T/I and t(k) = T/I are obtained from the itera- 
tive procedure. 

The equations (7) and (8) are nonlinear and require 
a self-consistent solution that is obtained using a conju- 
gate gradient algorithm. In this work, we are interested 
in analyzing all the stationary solutions for the current. 
As J(y)and V(J) are multivalued functions, the mul- 
tiple solutions for the quantities &o> b\, Ao, Aiare ob- 
tained by fixing, first the external potential V and af- 
terwards the current J in the equations of motion of the 
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system. This method permits us to obtain the complete 
J- ^characteristic curve. 

Once the amplitudes a^ a are known, the current is 
numerically obtained from, 

k,(7 

We study first a model which consists of two leads 
equally connected(Vt = Vr = V$) to two quantum 
dots with fi L = -V/2 and n R = V/2, t = 30 r , 
Vo = 5.48 r , £o = £i = —3.5 T (Kondo regime with 

t k w io~ 3 r with r = 27ry 2 /9 (o)). 

The Fig 1 shows the J-V characteristic curve for dif- 
ferent values of t c . Within a range of the external bias 
voltage values, there are three stationary solutions for 
the current. We can identify two regimes. One for 
r < t c < t c0 ,( t c0 is a constant to be analyzed below ), 
where the J — V characteristic curve has a Z-shape and 
another one for t c > t C Q where the curve becomes loop- 
shaped. An experiment that would consist in studying 
this system just by simply continuously modifying the ap- 
plied potential V would characterize the behavior of the 
current as bistable. At a critical value Vc^ there would 
be a discontinuous reduction of the current when V is 
increased and other Vc^(V~C| > Vc\) where the current 
would rise, (t c < t c0 ), or diminish , (t c > t c0 ), discon- 
tinuously when V is reduced. In this case the interior 
of the bistability remains inaccessible to the experiment. 
However, the internal shape could be reached, as it is for 
the case of double barrier heterostructrures, employing 
a load line in the experimental measurement that per- 
mits to go along all the points of the J — V characteristic 
curve. 

When t c > To, the Kondo levels £o and £iare splittcd 
into two molecular states, corresponding to the bond- 
ing and antibonding levels. They are located in, e± = 

{(e + £i)/2 ± y / (e -ei) 2 +4t2}/2. In the Fig. 2a and 
2b we plot function of V, for two values of t c . 

The main features are: a) For V = 0, the bonding and 
antibonding levels satisfies e± = ±t c . When the bias V is 
increased of the order of Tk , e± are kept almost constant 
and the dots reacts to the external potential in a coherent 
way. For higher values of V , £+ and £_ approach each 
other and go through the multistablc region. The right 
dot is connected to the left reservoir by the mediation of 
the left one through an effective interaction V^J that for 
energies in the vicinity of the Fermi level of the order of 
Tk, is approximately given by V^f^ — t c /2nVLp(0). The 
effective levels have a clear different behavior whether the 
right dot is interacting stronger with the left than with 
the right V^J* > Vr or the contrary, < Vr . For 

the case we are analyzing, (Vl — Vr — Vo), this can be 
controlled manipulating the value of t c . Defining i c oas 
the frontier between these two situation, in the case when 
t c < t c o, increasing the potential beyond the mutistablc 
region, the coherence between the dots is lost and e + and 



e_ converge towards their own chemical potential Hl(r)? 
as display in Fig. 2a. The interactions of the dots with the 
leads Vl and Vr are equally renormalized as can be seen 
in Fig 2c. The system possesses an splitted Kondo peak. 
One sub-peak with its weight essentially located at the 
left dot with energy at the neighborhood of [i^ and the 
other located at the right dot near \ir. For t c > t cn by 
the contrary, due to the stronger interaction between the 
dots, they maintain their coherence and result to be par- 
tially disconnected from the right lead. This is reflected 
in the behavior of the renormalized interaction of each 
dot with its respective lead as shown in Fig. 2d and in 
their localized levels e + and e_ , that both follow as 

V increases presented in Fig2b. The Kondo peak, almost 
pinned at \xl, is in this case equally distributed among 
the two dots. 

Next we will discuss the effect of the asymmetric site 
energies case, £o 7^ £\- We model the asymmetric site 
energies by setting £0 = (1 — «;)eo.£i = (1 + fc)eo, where 
k is the asymmetry parameter and eo = — 3.5 Tq . Figs 
3 shows the J-V characteristics for various values of k 
at t c = 1.3r and t c — 2.5r . It is seen that the region 
of multistability is very sensitive to the change of n and 
in particular it disappears for sufficiently large values of 
it. The coherent or incoherent behavior of the dots dis- 
cussed above can be controlled in this case changing the 
asymmetry parameter k. Fig 4 shows a map with the 
regions of multistability depending on the bias voltage 

V and k. It is seen that the width of the multistability 
region has lower and upper limits, which depend upon 
the asymmetry in the dot-lead coupling. Similar results 
were obtained for an asymmetric dot-lead coupling case, 
V L ± Vr. 

In this paper we have studied the J-V characteristic 
curve of a double-quantum-dot in the Kondo regime. It 
exhibits a multistable behavior that enhances as the in- 
terdot coupling tunneling increases. It becomes Z-shaped 
or loop-shaped depending on the various coupling tunnel- 
ing parameters, which control wether the dots respond 
coherently or incoherently to the external applied poten- 
tial. Moreover, we show that it is possible to manipulate 
the multistability in a sample by changing the asymmetry 
in the dot-lead coupling. 
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FIG. 1. J-V curve for a) t c = To, b)t c = 1.3T , 
c)t c = 1.6r and d) tc = 2.5r 



FIG. 2. Figures 2 a) and 2 b) level energies e± vs V at 
t c = 1.3Fo t c — l.GTo respectively. Figures 2 c) and d) the 
corresponding renormalization of the dot lead coupling 



FIG. 3. J-V curve for the asymmetric site energies case 
at a) t c = 1.3r with k=0.005 (solid line), 0.015 (dash 
line), 0.04 (dot line) and 0.06 (short-dash line), b)t c = 2.5F 
with k=0.04 (solid line), 0.08 (dash line), 0.12 (dot line) and 
0.16(short-dash line) 



FIG. 4. Map V versus k at a) t c = I.ITq, b)t c = 1.2r , 
c)t c = 1.6F ,d)i c = 2.5F 
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